Thermodynamics of Two Dimensional Magnetic Nanoparticles. 
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A two dimensional magnetic particle in the presence of an external magnetic field is studied. Equi- 
librium thermodynamical properties are derived by evaluating analytically the partition function. 
When the external field is applied perpendicular to the anisotropy axis the system exhibits a second 
order phase transition with order parameter being the magnetization parallel to the field. In this 
case the system is isomorph to a mechanical system consisting in a particle moving without friction 
in a circle rotating about its vertical diameter. Contrary to a paramagnetic particle, equilibrium 
magnetization shows a maximum at finite temperature. We also show that uniaxial anisotropy in 
a system of noninteracting particles can be missinterpreted as a ferromagnetic or antiferromagnetic 
coupling among the magnetic particles depending on the angle between anisotropy axis and magnetic 
field. 
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The remarkable development of experimental tech- 
niques grants increasingly more access to the hitherto 
mysterious nanometric world. In the particular case of 
magnetism, a recent experimental study has succeeded to 
investigate single clusters of around 3 nm diameter, using 
a special micro-SQUID device (I]]. Also, a wide variety 
of magnetic microscopies are under constant devel- 
opment to measure the magnetization process of single 
clusters in the nanometer range. In this context, the 
classical Stoner-Wohlfarth (SW) (5| model has been fre- 
quently revisited, to test wether or not such simplistic 
uniform rotation model really provides a realistic pic- 
ture of the magnetization reversal of a nanoparticlc |^|,^| . 
The SW model does not take explicitly into account ther- 
mal fluctuations, and therefore it is strictly valid for very 
low temperatures or very strong anisotropics On 
the other hand, when the anisotropies are rather small 
and/or the temperatures are sufficiently high, the so- 
called superparamagnetic approach must be employed. 
In such case, Neel |l(]] and Brown Jll| introduced the 
idea that the decay toward a thermal equilibrium state 
of the magnetic moment of a single-domain ferromagnetic 
nanoparticle is mediated by thermal fluctuations, follow- 
ing a simple Arrenhius switching probability QJl^l • Also 
in this case, further improvements of the model have been 
introduced year after year in order to include additional 
terms which are usually present in real samples, such as 
anisotropy 13- 15 a nd dipolar interaction among mag- 
netic nanograins [L6|~^9|. 

In this letter, we present a mechanical analog to a mag- 
netic moment in the presence of an external field perpen- 
dicular to its uniaxial anisotropy axis. This particular 
system exhibits a second order phase transition as a func- 
tion of the external field, where the order parameter is 
the magnetization. Its thermodynamics is also analyti- 



cally solved showing the behavior of the magnetization 
and susceptibility as a function of temperature. Finally 
we demonstrate that uniaxial anisotropy in an ensemble 
of non interacting SW particles can be missinterpreted 
as a ferro or antiferromagnetic coupling among the par- 
ticles. 

Let us consider a single magnetic domain particle with 
all its atomic moments rotating coherently on a plane in 
the presence of an external magnetic field. The particle 
is characterized by a constant absolute value of the mag- 
netization m = m s V, V being the volume of the par- 
ticle and m s its saturation magnetization. Well below 
the Curie temperature of the specific magnetic element, 
m s is independent of particle volume and temperature. 
The energy of this nanoparticle has three contributions: 
anisotropy energy, either due to shape, stress or crys- 
talline structure of the particle; Zeeman energy, owing 
to the presence of an external field, H, and the mag- 
netic dipolar interaction among the atoms composing the 
grain, and which will not be considered in the present 
work. Under these conditions the energy reads: 



E = -rh-H-KV 



(1) 



where k is the anisotropy constant and the unit vector 
n represents the easy magnetization direction. We will 
initially focus on the particular case in which the external 
field is perpendicular to the anisotropy axis. Then, the 
energy of the single particle can be written as: 
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E= — mH cos 9 mH a sin . 

2 

where we have defined the anisotropy field H a such that 
kV = ^mH a and represents the angle between magne- 
tization and external field. 
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A more suggestive form is obtained by making the en- 
ergy dimensionless as follows: 



£ = 



E 



= — cos fc> sin 

mH 2 H 



Now, let us consider a rigid circular ring of radius R 
rotating along its vertical diameter at angular frequency 
lu with a frictionless particle of mass m free to move on 
the ring, as illustrated in Fig. 1. The particle is under 
gravity, mg, and normal force N due to the ring, which 
is rotating, and therefore the mass describes a circular 
orbit due to the centripetal acceleration originated from 
the horizontal component of N. By choosing the zero 
level of gravitational energy at the center of the ring, the 
Lagrangian for this system reads 
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mR 2 9 2 + -muj 2 R 2 sin 2 + maR cos 9 




FIG. 1. Mechanical rotating system. 

Last equation can be rewritten as C — T — U , with T be- 
ing the kinetic energy and U an effective potential energy 
that includes the effects of gravity and the rotation of the 
system. Introducing a dimensionless form one finds that 
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Therefore the analogy between this mechanical system 
and the magnetic particle describes by Eq. 1 is clear. The 
centripetal acceleration, Rut 2 , plays the role of the mag- 
netic anisotropy whereas the gravity field is equivalent to 
the external magnetic field. Obviously, the downwards 
component can be further increased by means of the ap- 
plication of an additional force, and in such case the total 
force F divided by the mass of the particle represents the 
applied field. Notice that achieving the anisotropy field is 
equivalent to apply a downwards force which completely 
compensates the vertical component of the normal force 
on the particle. 

When we analyze this mechanical system, a second or- 
der phase transition appears. Analogously, analyzing the 
magnetization along the field as a function of the param- 
eter H a /H, a second order phase transition is clearly 



observed. In the latter case, the order parameter is the 
angle 9. Figure 2 shows the energy as a function of 9 for 
two different values of H a /H . Clearly, there is a smooth 
transition between one minimum in the energy landscape, 
for H a /H < 1 , to two minima, for H a /H > 1. The 
transition occurs for an external field strong enough to 
make the equilibrium state a magnetized one along the 
field. On the other hand, when the anisotropy field is 
stronger, there are two minima. They correspond to the 
states with non zero magnetization component along the 
uniaxial anisotropy axis together with the magnetization 
component along the external field. Then, there is a sec- 
ond order phase transition for H a /H = 1, in analogy 
with the well-known SW predictions. 
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FIG. 2. Reduced energy landscape for a magnetic sys- 
tem with anisotropy axis perpendicular to the external field. 
The energy is plotted against the angle 9 between the mag- 
netization of the system and the external field. We illustrate 
two representative regimes, as a function of the ratio H a /H . 
The broken curve illustrates results for H a / H — 0.5 and the 
continuous line for H a /H = 1.4. 
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FIG. 3. Equilibrium angle, 9, between magnetization and 
external field, as a function of the ratio H a /H. 



2 



Figure 3 illustrates the order parameter 9 as a func- 
tion of the ratio H a /H. The second order phase tran- 
sition is clearly evidenced. This figure also illustrates 
the position-frequency phase diagram for the mechanical 
analog depicted in figure 1. 

Now we investigate the thermodynamic of this particular 
system in the canonical ensemble. The partition function 
Z2D reads: 
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mH cos 9 + \mH a sin 2 6 
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This integral can be analytically solved and the result 
can be written as a series in Bessel functions of integer 
argument. 
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By making H a = we recover the ideal case of a param- 
agnetic system in 2D. Hereafter thermodynamics follows 
and it can be computed analytically. Magnetization and 
susceptibility of the particle along the field are respec- 
tively given by 
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These formula are evaluated using the recurrence rela- 
tions for the first derivative of the Bessel functions. 

Figure 4 illustrates the magnetization along the exter- 
nal field as a function of temperature for a system of non 
interacting particles with m = 1000/is each and exter- 
nal field H = Q.lkOe. At T = 29.6K the magnetization 
exhibits a maximum for H a /H — 30. This behavior is 
clearly understood by using the mechanical analog, where 
at fixed rotation speed u> and zero kinetic energy, the 
equilibrium angle 9 is larger than the value obtained by 
a small increasing of kinetic energy in the particle, due to 
the asymmetry of the energy landscape around a mini- 
mum. Further increasing in kinetic energy makes the an- 
gle increases up to an asymptotic final value of 6 = ir/2. 
In our magnetic system, larger kinetic energy corresponds 
to higher temperatures. Therefore, in equilibrium and for 
H a /H > 1, the magnetization exhibits a maximum at a 
temperature T such that kT = 0.152(H a /H - l)mH 

Figure 5 displays the magnetic susceptibility along the 
field direction as a function of temperature. The broken 
line depicts results for H a /H = 0.9 and the continuous 
line for H a /H = 30. In the region of the double well 
potential, H a /H > 1, one clearly observes a maximum 
in both magnetization and susceptibility curves. 
The maximum of the susceptibility for H a /H « 1 
approximately occurs at a temperature T max given by 
mH / kT max = 1.33, which corresponds to the value for a 
paramagnetic particle in 2D. For H a /H = 30, the maxi- 
mum occurs at mH/kT max = 1.0. 
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FIG. 4. Magnetization along the external field axis, in 
units of Hb , as a function of temperature for a system of non 
interacting particles with m = 1000/is each, and external field 
H = Q.lkOe. The continuous line, with magnetization scale 
at the left, illustrates results for H a = 3kOe. The broken 
line, with magnetization scale at the right, depicts results for 
H a = 0.09fcOe. 
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FIG. 5. Magnetic susceptibility along the external field, in 
units of /is/fcOe, as a function of temperature for a system of 
non interacting particles with m — IOOO^ib each, and external 
field H — Q.lkOe. The continuous line, with magnetization 
scale at the left, illustrates results for H a = 3kOe. The broken 
line, with magnetization scale at the right, depicts results for 
H a = 0.09fcOe. 

Finally, in order to explore the role of anisotropy, Figure 
6 illustrates the inverse susceptibility for different angles 
a between anisotropy axis and external magnetic field. 
For a parallel orientation between anisotropy axis and 
magnetic field we observe a susceptibility curve resem- 
bling a ferromagnetic system. 
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FIG. 6. Inverse Magnetic susceptibilities in units of H/m 
along the external field as a function of temperature for a sys- 
tem of non interacting particles with m = 1000/is each, exter- 
nal field H — O.lfcOe and H a = 3kOe. We show three curves 
depending on the angle a between anisotropy axis and exter- 
nal field as indicated inside the figure. The dashed straight 
lines are high temperature extrapolations. 

However, for perpendicular orientation between mag- 
netic field and anisotropy axis, we observe a typical 
Curie- Weiss antiferromagnetic-like behavior, where by 
extrapolating from higher temperatures, the zero cross- 
ing in the temperature axis occurs at negative values. 
When a = 45°, a simple Curie law is observed. Then, 
the extrapolated transition temperature changes contin- 
uously from a positive value, for parallel orientation, to 
a negative one, when the anisotropy axis is perpendic- 
ular to the external field. As a function of the relative 
angle between the anisotropy axis and external field, the 
net effect of the anisotropy field resembles a second order 
transition from a ferromagnetic to an antiferromagnetic 
coupling. 

In conclusion, a system of non-interacting magnetic 
particles with parallel anisotropy axes under the pres- 
ence of an external magnetic field applied perpendicular 
to the anisotropy axes exhibits (at T = OK) a second 
order phase transition as a function of the ratio between 
external field to the anisotropy field. The order parame- 
ter being the magnetization of the system. Such system 
is isomorph to a mechanical system consisting in a parti- 
cle moving without friction in a circle rotating about its 
vertical diameter. At finite temperature the thermody- 
namics of the system can be solved analytically and the 
magnetization along external magnetic field shows a max- 
imum at a finite temperature, contrary to that of a para- 
magnetic particle. The inverse susceptibility of a system 
of non interacting magnetic particles as a function of tem- 
perature demonstrates that uniaxial anisotropy acts sim- 
ilarly to effective ferromagnetic or antiferromagnetic cou- 
pling, depending on the angle between magnetic field and 



anisotropy axis. The presented results can be useful in 
the interpretation of magnetic data obtained in nanocrys- 
tallinc materials, specially novel artificially grown pat- 
terned media for high density magnetic recording, which, 
depending on the aspect ratio, can display a strong 2D 
character. 
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